An Asymptotically Stable Semi-Lagrangian scheme in the Quasi-neutral Limit
نویسندگان
چکیده
This paper deals with the numerical simulations of the Vlasov-Poisson equation using a phase space grid in the quasi-neutral regime. In this limit, explicit numerical schemes suffer from numerical constraints related to the small Debye length and large plasma frequency. Here, we propose a semi-Lagrangian scheme for the Vlasov-Poisson model in the quasi-neutral limit. The main ingredient relies on a reformulation of the Poisson equation derived in [4] which enables asymptotically stable simulations. This scheme has a comparable numerical cost per time step to that of an explicit Centre de Mathématiques Appliquées (UMR 7641), Ecole Polytechnique, F-91128 Palaiseau, France. [email protected] INRIA Nancy-Grand-EST (CALVI Project), and IRMA-Université de Strasbourg and CNRS, F-67084 Strasbourg cedex, France. [email protected] 1-Université de Toulouse; UPS, INSA, UT1, UTM ; Institut de Mathématiques de Toulouse; F-31062 Toulouse, France. 2-CNRS; Institut de Mathématiques de Toulouse UMR 5219; F-31062 Toulouse, France. [email protected] IRMA-Université de Strasbourg and CNRS, and INRIA NancyGrand-Est (CALVI Project), F-67084 Strasbourg cedex, France. [email protected]
منابع مشابه
Chapter 6 . Semi - Lagrangian Methods
Semi-Lagrangian (S-L for short) methods, also called quasi-Lagrangian, have been around since the early 1960’s. T.N. Krishnamurti is probably the first to use S-L method. The S-L scheme becomes really popular following Andre Robert’s work in the early 80’s (Robert 1982) when he combined the use of S-L with semi-implicit method (SLSI) which allows for large time step size because both schemes ar...
متن کاملDetermination of Stability Domains for Nonlinear Dynamical Systems Using the Weighted Residuals Method
Finding a suitable estimation of stability domain around stable equilibrium points is an important issue in the study of nonlinear dynamical systems. This paper intends to apply a set of analytical-numerical methods to estimate the region of attraction for autonomous nonlinear systems. In mechanical and structural engineering, autonomous systems could be found in large deformation problems or c...
متن کاملPositivity-preserving and asymptotic preserving method for 2D Keller-Segal equations
We propose a semi-discrete scheme for 2D Keller-Segel equations based on a symmetrization reformation, which is equivalent to the convex splitting method and is free of any nonlinear solver. We show that, this new scheme is stable as long as the initial condition does not exceed certain threshold, and it asymptotically preserves the quasi-static limit in the transient regime. Furthermore, we pr...
متن کاملAn Efficient Neurodynamic Scheme for Solving a Class of Nonconvex Nonlinear Optimization Problems
By p-power (or partial p-power) transformation, the Lagrangian function in nonconvex optimization problem becomes locally convex. In this paper, we present a neural network based on an NCP function for solving the nonconvex optimization problem. An important feature of this neural network is the one-to-one correspondence between its equilibria and KKT points of the nonconvex optimizatio...
متن کاملSolving a system of 2D Burgers' equations using Semi-Lagrangian finite difference schemes
In this paper, we aim to generalize semi-Lagrangian finite difference schemes for a system of two-dimensional (2D) Burgers' equations. Our scheme is not limited by the Courant-Friedrichs-Lewy (CFL) condition and therefore we can apply larger step size for the time variable. Proposed schemes can be implemented in parallel very well and in fact, it is a local one-dimensional (LOD) scheme which o...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
برای دانلود متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید
ثبت ناماگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید
ورودعنوان ژورنال:
- J. Sci. Comput.
دوره 41 شماره
صفحات -
تاریخ انتشار 2009